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Abstract. In the recent paper we calculated the 7-ray spectra from pulsar wind nebulae (PWNe), assuming that 
a significant amount of the pulsar rotational energy is converted into relativistic nuclei. These nuclei accelerate 
leptons which are responsible for most of the observed electromagnetic emission from PWNe. Small part of nuclei 
also interact with the matter of the supernova producing 7-rays, which can also contribute to the observed spectra 
of young nebulae. Here we calculate the spectra of neutrinos from the interaction of nuclei inside the nebula and 
the expected neutrino event rates in the 1 km 2 neutrino detector from: the Crab Nebula (PSR 0531+21), the 
Vela SNR (PSR 0833-45), G 343.1-2.3 (PSR 1706-44), MSH15-52 (PSR 1509-58), 3C58 (PSR J0205+6449), and 
CTB80 (PSR 1951+32). It is shown that only the Crab Nebula can produce the neutrino event rate above the 
sensitivity limit of the 1 km 2 neutrino detector, provided that nuclei take most of the rotational energy lost by 
the pulsar. The neutrino event rate expected from the Vela SNR is comparable to that from the Crab Nebula but 
these neutrinos are less energetic and emitted from a much larger region on the sky. Therefore it may be difficult 
to subtract the Vela SNR signal from the higher background of the atmospheric neutrinos. 

Key words, supernova remnants: pulsars: general - ISM: neutrinos: theory - radiation mechanisms: non-thermal 
- nebulae: Crab Nebula (PSR 0531+21); Vela SNR (PSR 0833-45); G 343.1-2.3 (PSR 1706-44); 3C58 (PSR 
J0205+6449); MSH15-52 (PSR 1509-58); CTB80 (PSR 1951+32) 



1. Introduction 

The possibility that hadronic processes can contribute to 
the observed 7-ray emission from the PWNe has been 
considered in recent years by e.g. Cheng et al. (1990) 
and Aharonian & Atoyan (1996). Following these works 
and earlier results (e.g. Berezinsky & Prilutsky 1978), 
we discussed a model for the production of radiation 
in hadronic processes with the application to the very 
young PWNe (Protheroe, Bednarek & Luo 1998, Bead 
& Bednarek 2002). In these papers, the 7-ray and/or neu- 
trino fluxes were estimated from the interaction of nuclei, 
injected by the pulsar, with the matter and soft radiation 
of the supernova remnant during its early phase of devel- 
opment, i.e. within the first few years. A general model 
like this has recently become very popular, finding its ap- 
plication in the estimation of neutrino fluxes from the 7- 
ray bursts in terms of the supranova model (e.g. Guetta 
& Granot 2002, Razzaque, Meszaros & Waxman 2002, 
Dcrmer & Atoyan 2003), which postulate the formation of 
a very fast pulsar during the supernova explosion (Vietri 
& Stella 1998). In another paper the contribution of nu- 
clei from the pulsar to the 7-ray (and neutrino) spectrum 
was also calculated for the best known example of the 
older supernova remnant, i.e. the Crab Nebula (Bednarek 
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& Protheroe 1997). Recently, the neutrino flux from the 
Crab Nebula has been independently estimated by Amato, 
Guetta & Blasi (2003). The authors predict the few to sev- 
eral neutrino events in a 1 km 2 detector per year from this 



Recently we have considered a more general model 
for the hadronic and leptonic processes inside the PWNe, 
which are based on the proposition that a significant part 
of the pulsar rotational energy is taken by relativistic nu- 
clei (Arons and collaborators, e.g. Arons 1998). Our pur- 
pose was to calculate the 7-ray spectra from the well- 
known PWNe. From the comparison of the calculations 
with the observations of these 7-ray nebulae, we derived 
some free parametrs of the considered model in order to 
predict the 7-ray and neutrino fluxes from other nebulae, 
which may become potential targets for the next gener- 
ation of 7-ray and neutrino telescopes. The calculations 
of the 7-ray spectra from leptonic and hadronic processes 
were reported in the accompanying paper by Bednarek 
& Bartosik (2003, BB03). The neutrino spectra are cal- 
culated in the present paper. Most of the earlier calcu- 
lations of the 7-ray fluxes from the PWNe, which refer 
mainly to the Crab Nebula (e.g. De Jager & Harding 1992, 
Aharonian & Atoyan 1995, De Jager et al. 1996a, Hillas 
et al. 1998), were done based only on the leptonic ori- 
gin for this emission. The calculations of the 7-ray fluxes 
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from other nebulae are available only for a few objects, 
e.g. nebulae around PSR 1706-44 (Aharonian, Atoyan & 
Kifune 1997), PSR 1509-58 (Du Plessis et al. 1995), Vela 
pulsar (De Jager et al. 1996b). 

The neutrino fluxes from a few PWNe have recently 
been re-scaled from the observed TeV 7-ray fluxes by 
Guetta & Amato (2002), based on the simple assump- 
tion that in all these nebulae the 7-rays above 2 TeV are 
produced in a decay of 7r°, which originate in hadronic 
processes. This assumption is in contradiction to the pop- 
ular opinion on the relative importance of the leptonic and 
hadronic processes in PWNe and also with the results pre- 
sented in this paper. 

2. The model for high energy processes in the 
pulsar wind nebula 

The details of the model for high energy processes in the 
PWNe are described in our first paper (BB03). For con- 
sistency we repeat here its main features. The evolution 
of the supernova remnant, containing energetic pulsar, is 
described following the main points of the picture con- 
sidered soon after the discovery of pulsars by Ostriker & 
Gunn (1971) and Rees & Gunn (1974). Let us assume 
that when the explosion occurs the expansion velocity of 
nebula at its inner radius is Vo,sn and its initial mass is 
-^o.SN- However this expension velocity can increase due 
to the additional supply of energy to the nebula by the 
pulsar and can also decrease due to the accumulation of 
the surrounding matter. We take these processes into ac- 
count when determining the radius of the nebula at the 
specific time, t, by using the energy conservation, 

where 

L cm {t) = B*R 6 s n 4 /6c 3 » 1.3 x 10 45 £ 12 p- s 4 erg s" 1 , (2) 

is the pulsar energy loss on emission of dipole electromag- 
netic radiation, f2 = 2ir/P, and the period of the pulsar 
P = 10~ 3 P ms s changes with time according to 

PL(t)=Pl ms + ^ xlO-W 12 , (3) 

where Po.ms is the initial period of the pulsar and B = 
10 12 -Bi2 G the strength of its surface magnetic field. The 
expanding nebula increases the mass from the surrounding 
medium according to 

4 , 

M SN (*) = M ,sn + -Tr Psur R^ ch (t), (4) 

where p S ur is the density of the surrounding medium and 
i?Neb is the outer radius of expanding envelope at the time, 
t, which depends on the expansion history of the nebula, 

i?Ncb = / V S x{t')dt'. (5) 

JO 

The expansion velocity of the nebula, Vsn(^)j and its den- 
sity of matter, pNcb = 3MsN(i)/47ri?jJj eb (i), at the time 



t have been found by solving the above set of Eqs. (1-5) 
numerically. 

The pulsar loses energy in the form of relativistic wind 
which extends up to the distance i? s h- At this distance, 
the pressure of the wind is balanced by the pressure of 
the expanding nebula. We estimate the location of this 
shock at the time, t, by comparing the wind energy flux, 
determined by L cm (Eq.[3J), with the pressure of the outer 
nebula, determined by the supply of energy to the neb- 
ula by the pulsar over the whole of its lifetime (Rees & 
Gunn 1974), 

^cm 

The assumption that the pulsar wind ram pressure is bal- 
anced by the total pressure inside the nebula is different 
from that made in the previous paper (BB03) in which 
the wind pressure is balanced only by the magnetic field 
pressure. In fact this modification changes the location 
of the pulsar wind shock for very young nebulae but do 
not have influence on the calculated neutrino fluxes from 
them since the escape conditions of heavy nuclei from the 
nebula do not change significantly. 

Knowing how magnetic field depends on the distance 
from the pulsar in the pulsar wind zone, we can estimate 
the strength of the magnetic field at the shock region from 

B sh = ^B pl J^) 3 ^, (7) 

where a is the ratio of the magnetic energy flux to the 
particle energy flux from the pulsar at the location of the 
pulsar wind shock, i? pu i and B pu \ are the radius and the 
surface magnetic field of the pulsar, respectively. The evo- 
lution of a with the parameters of the pulsar is found by 
interpolating between the values estimated for the Crab 
pulsar, ~ 0.003, and for the Vela pulsar, ~ 1 (see Eq. 16 
and below in Bednarek & Protheroe 2002 for details). 

It is likely that a significant part of the rotational en- 
ergy of the pulsar, transfered through the shock radius 
i? s h is in the form of relativistic heavy nuclei. In fact, such 
an assumption can explain morphological features of the 
Crab Nebula and also the appearance of extremely ener- 
getic leptons inside the nebula accelerated as a result of 
the rezonant scattering of positrons and electrons by heavy 
nuclei (Hoshino et al. 1992, Gallant & Arons 1994). From 
the normalization to the observations of the Crab pulsar, 
Arons and collaborators ( e.g. see Arons 1998) postulate 
that the Lorenz factors of iron nuclei in the pulsar's wind 
should be 

71 a 0.3Ze$ opcn /m iC 2 , (8) 

where m; and Ze are the mass and charge of the iron nu- 
clei, c is the velocity of light, and $ op0 n = \jE cra /c is 
the total electric potential drop across the open magneto- 
sphere. Eq. (8) postulates that the pulsar with a specific 
period and a surface magnetic field accelerates nuclei mo- 
noenergetically. As assumed in the Gallant & Arons (1994) 
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log (E / GeV) 

Fig. 1. The muon neutrino and antyneutrino spectra from 
the pulsar wind nebula at the distance of 2 kpc at a 
different time after the pulsar formation: t = 10 2 ,3 x 
10 2 ,10 3 ,3 x 10 3 , and 10 4 yrs. The initial period of the 
pulsar is 15 ms and the surface magnetic field 4 x 10 12 
G. The supernova has the mass of 3 M Q and the initial 
expansion velocity of 2000 km s _1 . The atmospheric neu- 
trino background within 1° of the source, ANB, is marked 
by the thin full curves, horizontal (upper) and vertical 
(lower) (from Lipari 1993). 




log (t / yr) 

Fig. 2. The muon neutrino event rate in 1 km detector 
from the pulsar wind nebula at the distance of 2 kpc as a 
function of the time after the pulsar formation for differ- 
ent initial parameters of the pulsar: initial period 15 ms 
and surface magnetic field B = 4x 10 12 G (full curve), 5 
ms and B = 4 x 10 12 G (dashed curve), and 15 ms and 
B = 2 x 10 13 G (dotted curve). The rates for the neutrino 
events from the nadir (absorption in the Earth) and the 
zenith directions are marked by the thin and thick curves, 
respectively. 



model for the Crab Nebula, the nuclei take a significant 
part, x, of the total rotational energy lost by the pulsar. 
The calculations of neutrino rates from the PWMe, pre- 
sented below, are done for the value of \ = 0.8. In the 
next section we integrate the injection spectra of nuclei 
over the activity period of the pulsar in order to obtain 
the equilibrium spectrum of the nuclei inside the nebula 
at a specific time after pulsar formation. 

3. Production of neutrinos by nuclei 

The nuclei injected into the nebula suffer adiabatic en- 
ergy losses, due to the expansion of the nebula, and 
rare collisions with the matter of the nebula. The most 
energetic nuclei can diffuse out of the nebula. We in- 
clude all these energy loss processes for the nuclei (adi- 
abatic losses, fragmentation, escape) following Bednarek 
& Protheroe (2002), and BB03, to calculate the equilib- 
rium spectra of different types of nuclei inside the nebula 
at the specific time after pulsar formation (see Sect. 4.1 
and Fig. 2 in BB03 for example equilibrium spectra of 
nuclei inside the nebula). 

The nuclei, with the equilibrium spectra, interact with 
the matter of the supernova remnant and produce neutri- 
nos via decay of pions. We estimate the neutrino spectra 
by simple re-scaling of the 7-ray spectra calculated for 
these hadrons, applying the scaling brake model proposed 
by Wdowczyk & Wolfendale (1987). The number of pro- 
duced neutrinos is on average two times larger than the 
number of 7-rays (the multiplicity of charged pions in re- 
spect to the neutral pions is larger by this factor) and 
the average energies of neutrinos are a factor of two lower 
than the 7-rays. As an example, we show in Fig. 1 the 



spectra of muon neutrinos (and antineutrinos) produced 
at a different times, 10 2 yr, 3 x 10 2 yr, 10 3 yr, 3 x 10 3 yr, 
and 10 4 yr after the formation of the pulsar at the dis- 
tance of 2 kpc. The initial parameters of the pulsar are: 
the period 15 ms and the surface magnetic field 4x 10 12 G, 
and the initial parameters of the nebula (supernova rem- 
nant) are: the mass 3M and the velocity 2000 km s _1 . It 
is assumed that supernova exploded in the medium with 
a typical density of 0.3 cm~ 3 . We also show the atmo- 
spheric neutrino background (ANB) expected within 1° 
of the source as calculated by Lipari (1993). As expected, 
the intensities of neutrino spectra strongly depend on the 
age of the nebula due to the lower densities of matter in- 
side the nebula. The spectra also shift to lower energies 
due to the adiabatic energy losses of nuclei, more efficient 
escape of higher energy nuclei from the nebula, and lower 
energies of the freshly injected nuclei inside the nebula by 
older pulsars. In Fig. 2 we show also the expected neu- 
trino event rate in a 1 km 2 detector from pulsars with 
different parameters. The predicted event rates are above 
1 per year if the age of the pulsar (and nebula) is shorter 
than ~ 10 3 - 3 x 10 3 yrs. 

4. Expected neutrino rates from specific nebulae 

We calculate the expected neutrino spectra at the Earth 
and the neutrino event rates in a 1 km 2 detector from the 
PWNe for which the 7-ray fluxes have been calculated in 
our previous paper (BB03). The initial parameters of the 
nebulae and the pulsars for these sources have been taken 
as applied in calculations of 7-ray fluxes. 

To estimate the neutrino rates we applied the muon 
neutrino detection probabilities by the 1 km 2 detector cal- 
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Fig. 3. Neutrino spectra (V M ) expected from different neb- 
ulae: Crab Nebula - PSR 0531+21 (thick full curve), 
Vela Nebula - PSR 0833-45 (thick dashed), G343.1-2.3 - 
PSR 1706-44 (thick dot-dashed), 3C58 - PSR J0205+6449 
(thick dot-dot-dot-dashed), MSH15-52 - PSR 1509-58 
(thin dotted) and with close high density medium with 
300 cnr 3 (thick dotted), CTB80 - PSR 1951+32 (thin 
dashed). In all cases, except PSR 1509-58, the parent su- 
pernova explodes in the medium with the density of 0.3 
cm~ 3 . The ANB is marked by the thin full curves. 

culated by Gaisser & Grillo (1987). The neutrinos arriving 
to the detector from the nadir direction can partially be 
absorbed by the Earth. This effect is significant for neu- 
trinos with energies > 100 TeV. We take it into account 
by applying the neutrino absorption coefficients calculated 
by Gandhi (2000). These absorption effects are not very 
important in the case of most of the considered objects. 

4.1. The Crab Nebula 

The observed parameters of the Crab Nebula and the pul- 
sar (the present radius ^2 — 3 pc, the expansion velocity 
2000 km s _1 , the period ~ 33.4 ms) can be fitted by as- 
suming that: (1) the pulsar was born with the initial pe- 
riod of 15 ms, the surface magnetic field 3.8 x 10 12 G, as de- 
rived from the observed period and period derivative with 
the assumption that the pulsar loses most of its energy on 
the dipole electromagnetic radiation; (2) the parent Crab 
Nebula supernova has 3M Q (consistent with the observed 
density of matter inside the nebula), and expends with 
the initial velocity of 2000 km s" 1 at the distance of 1830 
pc from the Crab Nebula (Davidson & Fesen 1985). For 
the applied parameters, we successfully described the ob- 
served 7-ray spectrum from the Crab Nebula (see BB03) 
as a self-consistent composition of the 7-ray spectra pro- 
duced by leptons (< 10 TeV) and hadrons (> 10 TeV). 

The muon neutrino and antyneutrino spectrum from 
the Crab Nebula, calculated for the above parameters, 
is above the atmospheric neutrino background within 1° 
of the source (see Fig. 3). However the neutrino rate in 
the 1 km 2 detector, about one event per year, is rela- 
tively low (see Table 1). This rate is similar to the neu- 
trino rate predicted in our earlier calculations (Bednarek 



& Protheroc 1997), which are based on somewhat deferent 
assumptions. In that model we assumed that neutrinos are 
produced in the interactions of nucleons with the matter 
of expanding supernova remnant. The nucleons were dis- 
solved from the heavy nuclei during their propagation in 
the inner pulsar magnetosphere. 

4.2. The Vela Supernova Remnant 

The Vela pulsar is inside the compact nonthermal nebula 
with the radius of ~ 7' and the extended Vela Supernova 
Remnant (Vela SNR) with the radius of - 3.5°. The ob- 
served parameters of the pulsar allow us to estimate its 
surface magnetic field, 4.5 x 10 12 G, and the characteristic 
age of 11.300 yrs. We model the Vela SNR assuming that 
the initial parameters of this object were similar to the 
Crab Nebula, i.e. the initial pulsar period of 15 ms, the 
mass of Vela SNR of 3 M Q and the expansion velocity of 
2000 km s _1 . The distance to the Vela pulsar is taken to 
be ~ 300 pc (Caraveo et al. 2001), although the older liter- 
ature suggested the value of 500 pc (e.g. Cha, Sembach & 
Danks 1999). For these parameters of the pulsar and neb- 
ula, we estimate its real age to be ~ 6000 yrs, consistent 
with the observed dimensions of the nebula. 

The muon neutrino spectrum, calculated from the Vela 
SNR, is just on the level of ANB within 1° around the 
source (Fig. 3). The expected neutrino event rate is sim- 
ilar to the Crab Nebula neutrino event rate (Table 1), 
but the neutrinos from the Vela SNR are less energetic. 
Moreover, due to the very large size of the Vela SNR on 
the sky, ~ 3.5°, the expected ANB for this source may 
be about an order of magnitude higher than marked in 
Fig. 3. Therefore the observation of a neutrino signal from 
the Vela SNR is unlikely. 

4.3. The Nebula around PSR 1706-44 

The pulsar PSR 1706-44, and its nebula (G 343.1-2.3), 
show close similarities to the Vela pulsar and SNR. It 
has also been reported at TeV 7-ray energies (Kifune et 
al. 1995, Chadwick et al. 1998). These pulsars have simi- 
lar characteristic ages and present periods. However, due 
to the lower surface magnetic field of the PSR1706-44, 
equal to 3.1 x 10 12 G, and a longer period, 102 ms, the 
real age of the pulsar has to be closer to its characteristic 
age 17.400 yrs. As in the previous modelling we assume 
that the initial period of this pulsar is also 15 ms and the 
expansion velocity of the nebula 2000 km s _1 . For these 
parameters of the pulsar and nebula, the consistency with 
the observed parameters of this object is reached for its 
real age 16.000 yrs. The distance to PSR 1706-44 is taken 
1.8 kpc (Taylor & Cordes 1993). 

The neutrino spectra calculated from PSR 1706-44 are 
about two orders of magnitude below the ANB (Fig. 3), 
and the predicted neutrino event rate in a 1 km 2 detec- 
tor, 4 x 10~ 3 per year, do not allow their detection in a 
reasonable time. These small values are due to the signifi- 
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Nebula 


Crab Nebula 


Vela nebula 


G 343.1-2.3 


3C58 


MSH15-52 


CTB80 


(pulsar) 


(PSR0531+21) 


(PSR0833-45) 


(PSR1706-44) 


(PSRJ0205+6449) 


(PSR1509-58) 


(PSR1951+32) 


neutrino rate 


1.3 (1.1) 


1.15 (1.05) 


4 x 10~ 3 


0.02 


3 x 10~ 3 


5 x 10~ 3 


(> 0.1 TeV) 










(1.0) 


(0.03) 



cantly larger distance and higher real age of PSR 1706-44 
in comparison to the Vela pulsar. 

4.4. The Nebula 3C58 around PSR J0205+6449 

The pulsar PSR J0205+6449 and its nebula 3C58 is some- 
times identified with the historic supernova in 1181 yr 
(Thorsett et al. 1992). However, in such a case, using the 
required initial period of the pulsar of <~ 60 ms and the ex- 
pension velocity of the nebula of <~ 5000 km s _1 (to fit the 
observed dimensions of the nebula), we were not able to 
describe successfully the observed electromagnetic spec- 
trum from this source in the broad energy range (BB03). 

Therefore, in order to fit the observed electromagnetic 
spectrum, we followed another proposition for the pa- 
rameters of this pulsar and nebula: the age of 5000 yrs 
(Bictenholz et al. 2001), which is consistent with the char- 
acteristic age (Murray et al. 2002), initial period of 15 ms, 
and expansion velocity of 1000 km s _1 (Fesen 1983). These 
values allow to fit the present parameters of this system. 
The surface magnetic field of this pulsar is estimated to 
be 3.6 x 10 12 G and its distance on 3.2 kpc (Roberts et 
al. 1993). However for these parameters, the neutrino spec- 
trum from PSR J0205+6449 is about two orders of magni- 
tude below the ANB (Fig. 3), and the predicted neutrino 
event rate (see Table 1) is not detectable. 

4.5. The Nebula MSH15-52 around PSR 1509-58 

The pulsar PSR 1509-58, with the characteristic age of ~ 
1700 yr, and its complex nebula MSH15-52, are sometimes 
identified with the historic supernova SN 185 (Thorsett 
1992). The distance to these objects is put in the range 
4.2 kpc (kinematics of H I, Caswell et al. 1975) up to 5.9 
kpc (dispersion measure, Taylor & Cordes 1993). We apply 
the value of 5.2 kpc. The pulsar has the present period of 
<~ 150 ms, and the surface magnetic field of 1.5 x 10 13 
G. Its present parametrs can be explained by assuming 
that the initial pulsar period was 15 ms and its real age 
is close to its characteristic age. However, in order to fit 
the observed dimensions of the nebula, 5' x 10', we have 
to assume that the expansion velocity of the bulk matter 
in this nebula is 5000 km s" 1 , which is much higher than 
that applied in the modelling of other nebulae. 

The neutrino spectrum calculated from this nebula is 
about two orders of magnitude below the ANB (Fig. 3) 
and the expected neutrino event rate in 1 km 2 detec- 
tor is very low (Table 1). However this source has been 
marginally detected by the CANGAROO telescope at TeV 



energies (Sako et al. 2000). We were not able to describe 
such a high level of the 7-ray emission from MSH15-52 
by only leptonic processes (BB03). However, the thermal, 
optical nebula RCW 89, containg H a line-emitting fila- 
ments with the density of ~ 5 x 10 3 cm" 3 , coincide with 
the NW component of the remnant MSH15-52 (Seward 
et al. 1983). Therefore, if the nuclei, injected by the pul- 
sar are captured by these high density filaments, then the 
reported level of 7-ray emission from this object can be 
explained by the 7-rays produced in interactions of nu- 
clei with the matter with average density of <~ 300 cm~ 3 . 
Then, the accompanying neutrino emission is comparable 
to the ANB (Fig. 3), and the expected neutrino event rate 
in a 1 km 2 detector is close to 1 event per year. 



4.6. The Nebula CTB80 around PSR 1951+32 

The pulsar PSR 1951+32 has a short period, 39.5 ms, a 
large characteristic age of 1.1 x 10 5 yr, and a relatively 
weak surface magnetic field of 4.9 x 10 11 G. It lies inside 
the supernova remnant CTB80, which consists of a 10' x 6' 
compact nebula and shell-like extended component with a 
diameter of 30'. The estimated dynamic age of the nebula 
9.6 x 10 4 yr (for the distance 2.5 kpc) matches the char- 
acteristic age of the pulsar (Koo et al. 1990). Therefore in 
our modelling, we assume that these objects arc at the dis- 
tance of 2.5 kpc and have the real age close to the pulsar 
characteristic age. The observed parameters of the pulsar 
and nebula are consistent with the initial period of the 
pulsar of 15 ms and the expansion velocity of the nebula 
2000 km s _1 . The neutrino spectrum calculated for the 
above parameters is also orders of magnitudes below the 
ANB (Fig. 3), and the predicted neutrino event rate is 
below the sensitivity of a 1 km 2 neutrino detector (Table 
1)- 

The nebula CTB80 is sometimes regarded as the rem- 
nant of the recent supernova in the year 1408, found in 
Chinese records (Strom et al. 1980, Wang & Seward 1984). 
If this is the case, then the initial period of the pulsar has 
to be very close to the one observed and the energy lost 
by the pulsar is very low. Assuming the expansion veloc- 
ity of the nebula qcual to 1000 km s _1 , to fit the observed 
dimension of the compact nebula, the expected neutrino 
event rate in this case is about an order of magnitude 
higher than estimated above (see the number in brackets 
in Table 1), but still below the sensitivity of the future 
neutrino detectors. 



6 



W. Bednarek: Neutrinos from the pulsar wind nebulae 



5. Conclusion 

According to the model discussed here for the high energy 
processes in PWNe, in which most of the observable 7-ray 
emission originates in leptonic processes and the accom- 
panying hadronic processes contribute only to the higher 
energy part of the spectrum, only young PWNe (with the 
age < 1 — 3 x 10 3 yr and at the typical distance of 2 kpc) 
should be detected by a 1 km 2 neutrino detector with 
the rate above 1 event per yr, i.e. above the atmospheric 
neutrino background within 1° of the source. Between all 
considered PWNe, which are observed in the TeV 7-rays, 
contains a young 7-ray pulsar, or are considered as likely 
objects for TeV 7-ray detection, only neutrinos from the 
Crab Nebula should be observed by the 1 km 2 with the 
event rate of 1 per yr. The neutrino event rate from the 
Vela SNR is comparable to that of the Crab nebula, but 
due to significantly lower energies of neutrinos and a much 
larger solid angle of the source on the sky, the signal from 
this object may not be extracted from the atmospheric 
neutrino background. The neutrino spectra from other 
considered nebulae (pulsars): G 343.1-2.3 (PSR 1706-44), 
MSH15-52 (PSR 1509-58), 3C58 (PSR J0205+6449), and 
CTB80 (PSR 1951+32), are significantly below the ANB, 
and the expected neutrino event rate is far below the sen- 
sitivity of a 1 km 2 neutrino detector in a reasonable time 
of observation. 

The neutrino rates, calculated in terms of the model for 
the high energy processes in PWNe considered here, are 
lower than those obtained by Guetta and Amato (2002), 
who rescaled the observed 7-ray fluxes above 2 TeV to 
the neutrino fluxes by simply assuming that all this 7- 
ray emission originates in hadronic processes. This as- 
sumption is in contradiction to the widely accepted lep- 
tonic model for the origin of the bulk emission from the 
WPNe. According to Guetta and Amato (2002), Crab , 
Vela, PSR1706-44, and PSR1509-58, should be easily de- 
tected by the 1 km 2 neutrino detector during one year. We 
have calculated the ratio of the neutrino luminosity to the 
7-ray luminosity above 1 TeV for the PWNe with reported 
TeV 7-ray emission. This ratio is equal to L v jL 1 ss 58% 
(the Crab Nebula), 22% (the Vela Nebula), and 0.5% (the 
nebula around PSR1707-44). Only in the case of the Crab 
Nebula the predicted neutrino luminosity should be close 
to the observed TeV 7-ray luminosity above 1 TeV. The 
neutrino fluxes can be comparable to the observed 7-ray 
fluxes even for relatively old nebulae if the high density 
medium is present close to the PWNa. As an example 
in Bednarek & Bartosik (2003), we consider the case of 
MSH15-52 (PSR 1509-58) which is close to the high den- 
sity medium (Seward et al. 1983). For this PWNa the ratio 
L v jL 1 is ~ 77%. 

In the case of all considered nebulae we assumed that 
the pulsars are born with the initial periods of 15 ms, as 
derived for the Crab pulsar. In fact, pulsars with shorter 
periods can also explain the main observed parameters 
of these nebulae. If the initial periods of these pulsars are 
significantly shorter then the expected neutrino event rate 



increases. As an example, in Fig. 2 we consider the case 
of the pulsar with the period of 5 ms which produces the 
neutrino event rate by a factor of 2-3 higher than predicted 
for the 15 ms pulsar, for the nebulae older than ~ 10 3 
yr. Therefore, significant reduction of the initial period of 
the pulsar does not influence the main conclusions of the 
paper concerning the detection of neutrinos from specific 
nebulae by a 1 km 2 detector. 

We assumed that the density of a medium in which the 
supernova exploded is 0.3 cm~ 3 . The pulsars surrounded 
by young supernova remnants, with similar parameters 
to the ones considered here, should produce much higher 
fluxes of neutrinos if they exploded in the high density 
medium. For example, the presence of the cloud with the 
average density of 300 particles cm -3 close to MSH15-52 
(PSR 1509-58) should allow detection of neutrinos and 
explain the marginal detection of TeV 7-rays from this 
object (Sako et al. 2000), as a result of radiation produced 
during interactions of hadrons with the matter. 
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